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Abstract 

Criterion for the Riemann hypothesis found by Baez-Duarte in- 
volves certain real coefficients defined as alternating binomial sums. 
These coefficients can be effectively investigated using Norlund-Rice's 
integrals. Their behavior exhibits characteristic trend, due to triv- 
ial zeros of zeta, and fading oscillations, due to complex zeros. This 
method enables to calculate numerical values of Ck for large values of 
k, at least to k = 4 ■ 10 s . 

We give explicit expressions both for the trend and for the oscilla- 
tions. The first tends to zero and is therefore, in view of the criterion, 
irrelevant for the Riemann hypothesis. The oscillations can be further 
decomposed into a series of harmonics with amplitudes diminishing 
quickly. Possible violation of the Riemann hypothesis would indicate 
that the amplitude of some high harmonic increases. 
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1 Introduction 



Several years ago a new expansion for the Riemann zeta function has been 
found by the author 0, 0: 

UJ ,-\2-> r(i-j) k'.-s-i^V 2> k k\ ( > 



k=0 



where 



_ r(fc + x) 
1 )k ~ r (x) 



and 



is the Pochhammer symbol (having such an unfortunate denotation and 
sometimes called the rising factorial). Pochhammer symbols are in fact poly- 
nomials in variable x with integer coefficients equal to the Stirling numbers 
of the first kind. Since tend to zero fast enough as k — > oo expansion (JTJ) 
converges uniformly on the whole complex plane. 

In fact, there exists a whole class of expansions similar to The crucial 
thing is to remove the single pole of £(s) in s = 1 which may be achieved 
either by multiplication by s — 1 (or by any other function having single simple 
zero at unity) or by subtraction of (or by any other function having single 
simple pole at unity). In the second case we get a series converging even 
faster: 



k=0 



where 



( fW+2 >-*Tr) (4) 

One has also the freedom of choosing the "node" points at which expres- 
sions (JTJ) or (J3J) give exact values, however the choice of even positive integers 
seems most natural, since A\~ may be expressed by Bernoulli numbers and 
appropriate powers of it avoiding values such as C(3)- 
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Finally, on can go all the way writing 



2(s-l) 



1+8 r (2T) " 1+s S-^ir * 



(5) 



with 



and 



A k :=J2(-iy( k )f(2j + 2) 



j=0 

'C'- 1 )^)- 1 + 1 - log(27T) 



s 

where f(s) is also regular on the whole complex plane. One can check that 
expansion (jHJ) gives exact values £(0) = — \ as well as C'(0) = — ^ log(27r) 
irrespective of how many terms in the series is taken into account. 

In 2003 Luis Baez-Duarte, investigating expansion (JTJ), found an interest- 
ing criterion for the Riemann hypothesis (RH) [2]. The crucial thing is to 
estimate the asymptotic behavior of certain real numbers defined as 

■ t: iy C),^:- (6) 

More precisely, RH is equivalent to the statement that 

c fc < k~ 3 / i+£ ,Ve > 0. (7) 

If 

c k « fc" 3 / 4 

then all zeros of £ are simple whereas regardless of the validity of the RH 
(unconditionally) 

c k « k^l\ 

As Baez-Duarte pointed out, criteria of this type, i.e. relating RH to values 
of zeta at integer points, were known for a long time, however this one is 
definitely simpler and well-fitted for numerical investigations. Such inves- 
tigations were performed numerically by the author in 2003 leading to the 
observation (which was surprising for us at that time) that the global behav- 
ior of coefficients c k may be split into a trend (which dominates for low k not 
exceeding roughly 10 4 ) and subtle oscillations superimposed on this trend. 
(Similar splitting has been found by the author in the case of Li criterion for 
the RH, however in this case the oscillations are much more chaotic [7j.) 
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2 Binomial transforms and asymptotics 



Of course, such behavior may be well understood in the theory of binomial 
transforms using Rice's integrals (already known by Norlund). Sums of this 
type can be interpreted as high order differences of appropriate numerical 
sequences {(p k }. They were carefully investigated see also [TQJ . The main 
result is (cf. Theorem 2 in j3j): 



k=riQ 



<p(k) 



■l) n ^Res 



<f(s) 



s fs-1 



is — n) 



(8) 



The proof of this fundamental theorem consists in applying classic Cauchy 
residue formula. 

In the case of © we have n = and (p(s) = l/£(2s + 2), which is analytic 
on [0, oo[ and meromorphic on the whole complex plane C, therefore it fulfils 
the assumptions of the above theorem. 



3 Results 

It is clear that the problem is to find all poles of suitable function, in our 
case l/£(s). These poles are due to simple zeros of zeta in s = —2n as well 
as due to complex zeros which we write as p = | (Obviously if RH were 
true then all 7s would be real.) It is a matter of elementary exercise to find 
that for n — 1, 2, ... the residues in real poles are: 

( 1 \ 1 (-l) n (2vr) 2n 

Res — -; s = -2n = = 2 / — 9 

\((sY J C'(-2n) (2n)!C(2n + l) KJ 

where the last equality is a consequence of the functional equation for the 
zeta function. 

Using main theorem (JHJ) and (JHJ) we find that the asymptotic form of the 
trend which stems from trivial zeros is 



c k 



k\ ^ (-l) m 7T 2m 1 . . 



71 



3/2 



m=2 
00 



T(k + m+ l)r(m - \) ((2m - 1 
B(k + l,m) (-l) m (27r) ; 



-T 

4tt 2 ^ r(2m- 1) C(2m- 1) 
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(the latter being numerically still more effective) whereas the oscillating part 
due to complex zeros in the critical strip is: 



V v 2 ; vc(2 S + 2)' 2 ; 

where B(x, y) is the Euler beta function and Ck = Ck + c. Both <fTTJ|> and (fTTf 
converge quickly, hence they are suitable for numerical estimations contrary 
to the direct approach using the main definition © which becomes very 
time-consuming as k grows. What's more, values of zeta function at positive 
even integers should be calculated with many significant digits. Since C{2n) 
tends quickly to unity as n grows it is advisable, when using JBj), to tabulate 
appropriate number of high-precision values of C(^ n ) — 1 to preserve suffi- 
cient amount of significant digits as well as to avoid repeated unnecessary 
calculations. Nevertheless, calculating of the single coefficient C300000 took 
about 2 weeks on a fast cluster of four computers ^2j- On the other hand, 
formulas (fTU|) and fTTf are much more effective: obtaining C1000000 is a matter 
of few tens of seconds on a modest machine. In numerical calculations using 
(II 1J) the function NResidue implemented in Mathematica proved especially 
useful. Figures 1-5 and the table below present the results. 



5 



k 




1 n 5 


j_i fif|Q7fi ■ 1 0~ 9 


2 • i n 5 


—7 89739 • 1 0~ 9 


3 • i n 5 


82876 • 1 0~ 9 


4 • i n 5 


—2 89364 • 1 0~ 9 


5 • i n 5 


—3 45567 • 1 f)~ 9 


6 • i o 5 


-4-1 1 3652 • 1 f)~ 9 


7 • i n 5 


-4-3 1 442Q • 1 f)~ 9 


8 • ID 5 


-1-2 00526 • 10~ 9 


9- 10 5 


-1.70316 • 10~ 10 


10 6 


-1.77502- 10~ 9 


2 ■ 10 6 


+8.08716- lO" 10 


3- 10 6 


-8.22419 • lO" 10 


4- 10 6 


+8.01923- 10~ 10 


5- 10 6 


+2.78245 • 10~ 10 


6- 10 6 


-5.00102 • 10~ 10 


7- 10 6 


-5.21564- 10~ 10 



The residues in (|llj) may be expressed by complex zeros of ( in a manifest 
way. Introduce function £(s) = £(1 — s) as usual 

((s) := 2(3 - l)n-^T (l + |) C(«). 

Since £ may be factorized using Hadamard product 

<M-nH) 

(the product being taken over paired complex zeros p) one can see that the 
residue of l/£(s) in each particular pole p may be expressed as a function of 
all remaining ps. Let us label consecutive complex zeros as p i according to 
the increase of |7|. Choose further a particular zero p n and introduce certain 
"crippled" function £ fc as 

£(«) = (s-p n )£ n (s) (12) 

: — i( 1 -s)n(i-a # n(i-a 
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l/£ fe (s) is obviously regular at p k and 




Res ( 1 ■ 8 = -I±M = - 
lenotes complex conjugation and 



2£n(Pn) 



1 h ( P * n ) 



(13) 



/i(s) :=2(s-l)7r- s/2 r + | 



) 



which is regular at p. Formula (|T3~j) is not very useful in practice since the 
product over paired complex zeros converges slowly, nevertheless, when com- 
bined with (fTTj). gives an explicit formula relating any Ck to complex zeros 
P- 

Using different approach Baez-Duarte gave another version of explicit 
formula of this kind, cf. [2], Theorem 1.5. The present approach is based 
entirely on Norlund-Rice's integral and provides, I believe, immediate and 
more natural interpretation for the existence of oscillatory component of the 
coefficients c k which is completely hidden in their primary definition (jUJ). 
After all, the behavior of this very component is crucial for the RH, see Fig. 
5. Last not least, resulting formulas are more handy and more effective in 
numerical calculations. 

4 Discussion and open questions 

Having any criterion for the RH the key thing is to say how "useful" it may 
be in numerical experiments. For example, Li's criterion states that RH is 
true if certain numbers A n are positive [7j. However, the fact that n initial A n 
are positive implies that roughly only y/n complex zeros lay on the critical 
line . Since high A n are extremely difficult to compute we can honestly say 
that Li's criterion, although very elegant, is pretty useless in practice. The 
natural question arises now: does the fact that k initial c k obeys (J7J) implies 
that certain number m of initial complex zeros lay on the critical line? What 
is the relation between k and ml Below we present simple argument based 
on formulas derived in this paper. 

The consecutive residues (fT3~]) seem to be bounded sequence of complex 
numbers, cf. Fig. 6, so we neglect their influence. Using Euler beta function 
expansion for large values of its first argument ^H] 



B(a,b) oc T(b)a- b 1 



b{b~l) 



2a 



( 




) 



) 



a 



(14) 
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we get: 



B(k + i, 1 -±p)oc(k + i)-^r( 1 -±l)= (15) 



= r 0-^) (* + i)" 1 [ cos (I lo S ( k + !)) + i sin (| log (k + 1))] (16) 

It is clear that if 7 is real then we get bounded oscillations in and, in the 
view of (JTj), RH is satisfied. Suppose, however, that there exists somewhere 
defiant extremely high zero | + ij with 7 having non- vanishing imaginary 
part. (More precisely, there should be four such zeros at once placed sym- 
metrically with respect to the critical axis and the real axis.) We are now 
sure that such a zero, if any, lies higher than 7 ~ 10 13 [5], ^1] and probably 
higher than 7 ~ 10 21 jSj. This non- vanishing imaginary part would cause 
the trigonometric functions in (|16[) to acquire, roughly speaking, growing 
amplitudes thus violating criterion (JJJ). Now, the amplitude of such grow- 
ing oscillations would be extremely small due to elementary properties of 
the Euler beta function which appears in (|llj). Specifically, for 7 = 2- 10 21 
(which is not much higher than the range acquired recently by Odlyzko) the 
amplitude due to the corresponding zero is roughly io -3 - 510 . Of course, 
the task to extract this particular harmonic, let alone to check whether its 
amplitude grows, is far beyond any numerical capabilities. In other words, 
even if the relatively low complex zero of zeta, say the hundredth, had some 
small shift off the critical line we would not be capable of finding this using 
© since its amplitude, in the sense of (fTBj). is about 10 78 times smaller than 
the amplitude of the first zero. 



This is rather sad news for all those who wish to (dis)prove the cele- 
brated RH using simple criteria. On the other hand, estimations like that 
presented above enable one to understand why this problem is so tough. The 
Holy Grail and immortality it would bring to persistent searcher remain still 
out of reach... Well, perhaps there exists even simpler explanation to this 
mystery, explanation not involving all that terrible mathematics. All things 
work together for good to them that love God says the inscription on Bern- 
hard Riemann's tomb taken from St. Paul's Epistle to the Romans (8:28). 
Shouldn't these words be taken for a clue? Yet, for a century and a half 
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many mathematicians would sell their souls for the proof. Clearly, there is 
an apparent contradiction between such a deal and St. Paul's words. 

Figure captions 

Fig. 1. Coefficients Ck for k up to 10 4 exhibit no apparent sign of oscil- 
lations. After calculating numerically the first 10 3 coefficients Ck this looked 
clS ct pleasant smooth curve" 

Fig. 2. Million coefficients Ck computed using (fTU|) and (JTTJ). Only 6 
initial terms of the series were needed to accomplish sufficient accuracy. 

Fig. 3. The same as Fig. 2 with logarithmic scale in k. 

Fig. 4. Behavior of Ck (yellow) may be split into strictly growing trend 
(red) plus the oscillating part (blue) here plotted in the logarithmic £;-scale. 
Initially the trend dominates over the oscillations but it tends to zero faster 
than these oscillations. 

Fig. 5. Over 44 million components of Ck- trend Ck (red) and oscillating 
part Ck (blue), both multiplied by fc 3 / 4 (compare RH criterion (J7J)). The 
former grows strictly to zero whereas the latter seems to tend to perfect 
sine wave which is dominated by the first zero. (Amplitudes of further zeros 
diminish quickly.) If this tendency persisted to infinity for all zeros RH 
would be true. Maximum value of k — 44700000 is related to parameter 
$MinNumber in a particular version of Mathematica. 

Fig. 6. 600 initial residues (JT3J). Red color denotes real part, green - 
imaginary part. Points are joined together in order to better visualize their 
behavior. 

References 

[1] L. Baez-Duarte, On Maslanka's Representation for the Riemann Zeta 
Function, preprint, 11 July 2003; this preprint fills the gap in an original 
proof given in (0). 

[2] L. Baez-Duarte, A Sequential Riesz-Like Criterion for the Riemann Hy- 
pothesis, International Journal of Mathematics and Mathematical Sci- 
ences 2005:21 (2005) 3527-3537. 

[3] L. Baez-Duarte, A New Necessary and Sufficient Condition for the Rie- 
mann Hypothesis, preprint, July 13, 2003 (earlier version of [2]). 



9 



[4] Ph. Flajolet, R. Sedgewick, Mellin Transforms and Asymptotics: Fi- 
nite Differences and Rice 's Integrals, Theoretical Computer Science 144 
(1995), no. 1-2, 101-124. 

[5] X. Gourdon, The 10 13 first zeros of the Riemann Zeta function, and 
zeros computation at very large height, preprint, October 2A th 2004. 

[6] K. Maslanka, The Beauty of Nothingness: Essay on the Zeta 
Function of Riemann, Acta Cosmologica XXIII-1, p. 13-18, 1998; 
A hypergeometric-like Representation of Zeta function of Rie- 
mann, posted at arXiv:math-ph/0105007 vl 4 May 2001; see also 
http://functions.wolfram.com, citation index: 10.01.06.0012.01 and 
10.01.17.0003.01. 

[7] K. Maslanka, Li's Criterion for the Riemann Hypothesis - Numerical 
Approach, Opuscula Mathematica, vol. 24 (2004), no. 1, s. 103-114. 

[8] A. M. Odlyzko, The 10 22 -nd Zero of the Riemann Zeta Function, Con- 
temporary Mathematics. 

[9] J. Oesterle, Paris, uncirculated note, circa 2001 (information from A. 
Voros) . 

[10] Th. Preu, Rice's integrals in: 2nd Joint Advanced Student School 2004, 
Complexity Analysis of String Algorithms, St. Petersburg, March 28 - 
April 7, 2004 (pedagogical lecture based on jl]). 

[11] S. Wedeniwski, ZetaGrid - Verification of the Riemann Hypothesis, 
available at http://www.zetagrid.net/. (The project started in August 
2001 and was closed in December 2005.) 

[12] M. Wolf, private communication, 2006. 

[13] http://functions.wolfram.com/06.19.05.0007 



10 



Fig. 1 
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